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Abstract Fieldwork in Jiemuxi National Nature Reserve, Hunan, China, discovered morphologically indis- 
tinguishable groups of brown frogs that bred at different, exclusive times of the year. A principal components analysis 
of morphometric data, molecular analyses, and the exclusive breeding season suggested the occurrence of two species. 
The population that breeds during the winter was found to be an undescribed species and was subsequently described 
as Rana jiemuxiensis sp. nov. The new cryptic species can be identified from its congeners at the same locality by 
having a different breeding season and a divergent DNA barcode. Six major lineages of brown frogs were resolved 
with high nodal support. Rana japonica, R. chaochiaoensis, R. omeimontis, and R. jiemuxiensis have independent 
maternal lineages. Rana hanlunica and R. maoershanensis have essentially identical maternal lineages and they appear 
to represent the same species. The validity of the species R. Jongicrus, R. zhenhaiensis, and R. culaiensis and their 
genealogical relationships are not resolved and deserve further study. The genealogy indicates that sympatric cryptic 
lineages probably widely exist in the R. Jongicrus group. This highlights the importance of future more fine-scaled 
samplings and the inclusion of nuclear genes. 
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similarities make many species notoriously difficult 
to identify (Che et al., 2007; Liu and Hu, 1961), 
especially when closely related species have overlapping 


1. Introduction 


Brown frogs of the genus Rana Linnaeus, 1758 (type 
species: Rana temporaria Linnaeus, 1758), also known 
as wood frogs, are mainly distributed from temperate 
Eurasia into Indochina, and western North America (Frost, 
2010). This genus contains 48 species. Morphological 


distributions (personal observation). 

Presently, 22 species of Rana are recognized from 
China (Frost, 2010; Yang et al., 2010), six of which 
are assigned to the R. Jongicrus group by Fei et al., 


(2009) (also called the R. japonica group) as follows: 
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hanluica Shen, Jiang et Yang; R. longicrus Stejneger; 
R. omeimontis Ye et Fei, and R. zhenhaiensis Ye, Fei 
et Matsui. In addition, two recently described species, 
R. culaiensis Li, Lu et Li and R. maoershanensis Lu, 
Li et Jiang, are considered to be junior synonyms of R. 
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zhenhaiensis and R. chaochiaoensis, respectively (Fei et 
al., 2009). 

The R. Jongicrus group is widely distributed in 
southern China including Taiwan. Pope and Boring (1940) 
consider most of the species as R. japonica Boulenger 
1879. The recent surge in new species descriptions (Li et 
al., 2008; Lu and Li, 2002; Lu et al., 2007; Shen et al., 
2007) suggests that diversity within the group may be 
underestimated, especially in poorly investigated forested 
regions of southern China. The morphology appears to 
be highly conserved, assuming that the group contains 
many cryptic species. Alternatively, many of the newly 
described species may not be valid. The reconstruction 
of their evolutionary history should shed light on the 
problem. Herein, DNA barcoding with the cytochrome 
c-oxidase subunit I (COI) (Hebert et al., 2003) and 
cytochrome b (Cyt b) genes is used to identify genetically 
distinct lineages worthy of more intense taxonomic 
study (Smith et al., 2008) and to hypothesize the groups’ 
matrilineal history. 

During recent herpetofaunal surveys in Hunan 
province, we collected brown frogs in Jiemuxi National 
Nature Reserve during two mutually exclusive breeding 
seasons at one locality. This suggested either the 
presence of two breeding seasons for the species, or 
two cryptic species with different breeding seasons. The 
frogs that bred at different seasons could not be readily 
distinguished morphologically. To differentiate between 


these possibilities, we formed a testable hypothesis. Our 
null hypothesis of conspecificity could not be rejected 
by DNA barcodes alone because mitochondrial DNA 
cannot document the absence or presence of gene flow. 
The null hypothesis could be rejected by using DNA 
barcodes as an indicator of distinctiveness, matrilineal 
history to infer the co-occurrence of lineages in breeding 
frogs, and morphological and ecological data as indicators 
of niche partitioning. For example, the presence of 
exclusive lineages during a particular breeding season 
would indicate lineage segregation and thus serve to 
reject the null hypothesis. Either the absence of discrete 
lineages or the sharing of lineages in the two breeding 
seasons would fail to reject the null hypothesis. Similarly, 
morphological distinctiveness can also be used to reject 
the null hypothesis. The investigation involved all the 
other species in the R. longicrus group widely distributed 
in southern China, except for R. chevronata, and Japan. 


2. Materials and Methods 


2.1 Sampling From 2008 to 2010, field surveys were 
conducted in southern China covering the entire range 
of the R. Jongicrus group including Shandong, Henan, 
Jiangsu, Zhejiang, Anhui, Jiangxi, Hubei, Hunan, Fujian, 
Guangdong, Guangxi, Yunnan, Sichuan, and Guizhou 
provinces plus Chongqing Municipality. Adults, juveniles, 
and tadpoles were collected. Following euthanization, 


Figure 1 Geographic distribution of the ingroup taxa from southern China used in this study. See Table 1 for information on species and 


sampling sites. 
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Table 1 Voucher specimens, localities, and GenBank accession numbers for Chinese brown frogs, Rana. GBN = GenBank Accession 


No. KIZ = Kunming Institute Zoology, CAS. SCUM = Museum of Sichuan University. * = type locality. 


Species Voucher No. Locality No. Locality COI Cyt b 
GBN GBN 
Ingroup 
R. chaochiaoensis SCUM045096WD 1 Zhongdian, Yunnan, China JF939066 JF939106 
KIZ06425 2 Zhaojue, Sichuan, China* JF939103 JF939140 
R. culaiensis KIZSD080501 12 Mt. Culaishan, Shandong, China* JF939082 JF939119 
KIZSD080502 12 Mt. Culaishan, Shandong, China* JF939083 JF939120 
KIZ03407 7 Mt. Maoer'shan, Guangxi, China JF939070 JF939110 
K1Z03408 7 Mt. Maoer'shan, Guangxi, China JF939071 JF939111 
KIZYPX10630 7 Mt. Maoer'shan, Guangxi, China JF939085 JF939122 
KIZYP0809001 4 Lichuan, Hubei, China JF939084 JF939121 
KIZYPX1142 8 Mt. Yangmingshan, Hunan, China JF939095 JF939132 
KIZYPX1144 8 Mt. Yangmingshan, Hunan, China JF939096 JF939133 
KIZYPX10775 6 Xinhua, Hunan, China JF939088 JF939125 
KIZYPX10776 6 Xinhua, Hunan, China JF939089 JF939126 
KIZ02409 11 Mt. Jigongshan, Henan, China JF939068 JF939108 
R. jiemuxiensis K1IZ05263 5 Jiemuxi National Nature Reserve, Hunan, China* JF939090 JF939127 
KIZ05264 5 Jiemuxi National Nature Reserve, Hunan, China* JF939091 JF939128 
KIZYPX10817 5 Jiemuxi National Nature Reserve, Hunan, China* JF939092 JF939129 
R. hanluica KIZ03477 9 Mt. Mangshan, Hunan, China JF939098 JF939135 
KIZYPX1172 8 Mt. Yangmingshan, Hunan, China* JF939099 JF939136 
KIZYPX1177 8 Mt. Yangmingshan, Hunan, China* JF939100 JF939137 
KIZYPX1182 8 Mt. Yangmingshan, Hunan, China* JF939102 JF939139 
KIZHuN024 5 Jiemuxi National Nature Reserve, Hunan, China JF939076 JF939114 
KIZHuNO1 5 Jiemuxi National Nature Reserve, Hunan, China JF939075 JF939113 
KIZHuN06 3 Jiemuxi National Nature Reserve, Hunan, China JF939077 JF939115 
KIZGX84 7 Mt. Maoer'shan, Guangxi, China JF939074 JF939112 
KIZYPX10654 7 Mt. Maoer'shan, Guangxi, China JF939086 JF939123 
KIZYPX10656 7 Mt. Maoer'shan, Guangxi, China JF939087 JF939124 
R. longicrus KIZ15026 16 Nanzhuang, Miaoli, Taiwan, China JF939067 JF939107 
KIZYPX11133 10 Shaoguan, Guangdong, China JF939093 JF939130 
KIZYPX11138 10 Shaoguan, Guangdong, China JF939094 JF939131 
KIZYPX11472 15 Fuzhou, Fujian, China JF939097 JF939134 
KIZYPX14268 14 Fuding, Fujian, China JF939104 JF939141 
R. omeimontis KIZ02424 3 Mt. Emeishan, Sichuan, China* JF939069 JF939109 
KIZSC148 3 Zhangceun, Sichuan, China JF939081 JF939118 
R. zhenhaiensis — K1Z0803271 13 Zhenhai, Zhejiang, China* JF939065 JF939105 
R. japonica KIZYPX11775 Japan JF939101 JF939138 
Outgroup 
R. amurensis KIZ070423558 Shangzhi, Heilongjiang, China JF939079 JF939117 
R. chensinensis | KIZRDOSSHX01 Huxian, Shaanxi, China JF939080 Zhou et al. unpubl. 
data 
R. dybowskii KIZ070423448 Huangnihe, Jilin, China JF939078 JF939116 
R. huanrensis SYNU040006 Huanren, Liaoning, China JF939072 Zhou et al. unpubl. 
data 
R. kukunoris CJ06102001 Qinghai Lake, Qinghai, China JF939073 Zhou et al. unpubl. 
data 


tissues dissected from adult specimens were preserved 
in 95% ethanol. Voucher specimens, including juveniles 
and tadpoles, were fixed in 10% buffered formalin, and 
all metamorphosed frogs were later transferred to 70% 
ethanol. All specimens were deposited in the Kunming 
Institute of Zoology, Chinese Academy of Sciences 
(KIZ). Comparative tissues from Taiwanese samples were 
provided by the National Museum of Natural Science 
(NMNS) Taiwan, China. Sampling data, which involved 


type localities, were summarized in Table 1 and Figure 
1. Rana chensinensis, R. huanrenensis, R. kukunoris, R. 
dybowskii, and R. amurensis were used as outgroup taxa in 
the genealogical analyses according to their phylogenetic 
positions as hypothesized by Che et al. (2007). 


2.2 Morphometrics A total of 109 specimens including 
five species were examined (Table 2); R. chevronata has 
not been seen since its description. Measurements were 
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Table 2 Linear measurements (in mm) of five Chinese species of brown frogs, genus Rana used in principal components analysis. The 


feature abbreviations are provided in text 


R. chaochiaoensis R. hanluica 


R. jiemuxiensis 


R. omeimontis R. zhenhaiensis 


Feature Range Mean+SE Range Mean+SE Range Mean+SE Range Mean+SE Range Mean+SE 
SVL 40.22-59.91 49.89+6.94 46.41-68.58 58.9747.66 34.11-53.57 44.00+4.91 49.37-62.60 56.2743.20 30.16-51.63 41.2747.21 
HL 12.12-18.51 15.84+2.09 14.61-25.16 18.54+2.65 10.37-17.38 13.3141.74 16.37-19.48 17.57+0.85 —9.34-16.52 12.88+1.99 
HW 11.97-21.22 16.5042.74 14.47-22.26 18.62+2.40 10.02-16.80  13.9141.58 15.38-20.03 17.76+1.41 9.27—16.28 13.08+2.17 
SL 6.64-9.36 7.97+0.88 7.47-11.27 9.501.21 4.05-12.50 7.0141.30 7.57-9.29 8.47+0.56 4.28-8.16 6.66+1.23 
EL 4.37-6.37 5.1140.59  4.26-6.57 —5.7140.60  2.82-6.16 4.50+0.64 4.92-6.92 5.7240.50 2.84-5.22 4.39+0.77 
IND 2.72-4.70 3.6740.54 3.40-5.21  4.2540.42 2.54-3.81 3.2040.33  3.42-4.71 4.00+0.33 2.59-3.98 — 3.20+0.43 
IOD 5.16-8.24 6.67+0.95 6.44-9.67  7.85+0.84 4.31-7.50 5.74+0.74 5.91-7.76 7.06+0.42 4.31-6.85 5.48+0.84 
UEW 2.96-5.24 3.7340.56 3.15-5.61 4.50+0.69 2.54-4.08 3.1940.41 3.28—4.73 4.19+0.40 2.68—4.62  3.56+0.55 
TYE 3.00-5.53 4.45+0.85 3.18-5.94  4.8340.64 2.48-4.11 3.45+0.47 3.98-5.59 4.92+0.40 2.51—4.13  3.46+0.45 
LAL 18.84-26.04 22.66+2.61 20.89-31.09 26.35+3.15 15.45-27.25 21.67+2.61 23.57-28.52 26.35+1.36 13.24-22.66 19.11+2.76 
FLI 5.08-7.68 6.18+0.84 5.27-9.19  7.09+1.15 3.78-7.26 5.1740.84 4.53-8.35 6.15+0.88 3.35-6.62 5.29+1.14 
FL2 4.12-6.90 5.2940.75 4.48-7.53  5.93+0.90 2.77-6.76 4.52+0.82 4.40-7.25 5.2840.69 3.10-5.70  4.61+0.80 
FL3 6.45-9.57 7.53+0.82 6.93-12.04 8.94+1.16 5.40-9.37 6.90+0.94 6.83-9.40 8.1740.80 5.40-8.93 7.31+1.08 
FL4 4.34-6.65 5.3840.73 4.33-7.25 5.83+0.79 3.16-6.32 4.46+0.63 4.35-6.78 5.5740.70 3.26-5.87  4.55+0.74 
FEL 22.17-35.22 29.91+4.72 27.78—43.75 36.51+4.52 18.27-31.51 25.48+3.39 28.57-37.81 33.98+2.40 15.20-30.44 22.93+4.51 
TL 20.29-37.29 31.39+45.11 10.70-47.26 38.43+7.29 19.98-34.77 27.73+3.82 31.18-40.55 36.60+2.44 17.40-31.96 24.90+4.59 
FTL 12.86-20.04 16.8142.36 16.36-22.98 19.94+1.96 11.37-18.92 14.78+1.89 17.25-21.74 19.50+1.17 9.71—-16.49 13.63+2.14 
FOL 22.99-34.07 28.47+3.46 27.58-36.87 33.06+2.71 15.54-34.58 24.78+3.99 30.85-35.72 33.05+1.60 17.32-30.62 24.71+3.91 
TEL 2.38-3.46 2.6740.28 2.18—4.39  3.5240.56 1.75-3.51 2.69+40.47  2.99-4.44 3.63+40.43 1.72-3.35 2.36+0.56 
TL2 3.11-6.56 4.5840.87 4.20-7.20 = 5.7140.79 2.96-7.12 4.4640.85 4.98-6.24 5.6640.40  2.40-4.51 3.87+0.70 
TL3 6.88-10.77 8.9441.18 8.23—-12.76 = 10.8+1.28 3.12-10.13 = 8.06+1.47  9.69-13.36 11.03+0.92 5.39-9.53 7.71+1.23 
TL4 12.79-19.95 16.24+2.21 15.43-24.00 19.4742.17 10.30-19.58 15.22+2.11 18.16-22.40 20.04+1.20 9.22-17.23 14.34+2.55 
TL5 7.58-12.73 9.58+1.47 9.01-14.36 12.08+1.56 4.12-11.66 8.741.51 10.47-14.38 12.11+0.95 5.98-10.03 8.21+1.17 


made by one person using dial calipers with a precision 
of 0.01 mm. Twenty-three linear measurements were 
taken as follows: SVL (snout-vent length), HL (head 
length), HW (head width), SL (snout length), EL (eye 
length), IND (internarial distance), IOD (interorbital 
distance), UEW (upper eyelid width), TYE (tympanum 
outer diameter), LAL (lower-arm length), FL1 (first finger 
length), FL2 (second finger length), FL3 (third finger 
length), FL4 (fourth finger length), FEL (femur length), 
TL (tibia length), FTL (length of foot and tarsus), FOL 
(foot length), TL1 (first toe length), TL2 (second toe 
length), TL3 (third toe length), TL4 (fourth toe length), 
and TL5 (fifth toe length). Because only three of the 23 
morphological characters differed significantly between 
females and males (TYE, FL2, and FL3; t-test, P < 
0.05), data from females and males were combined for 
a principal components analysis (PCA) implemented in 
SPSS 13.0. 


2.3 Extraction, amplification and sequencing Total 
genomic DNA was extracted from liver or muscle 


tissue using the standard phenol—chloroform protocols 
(Sambrook et al., 1989). A fragment of mitochondrial 
cytochrome oxidase subunit I (COI), the standard 
barcoding marker, was used to amplify all samples, 
together with Cyt b. Amplification was performed in a 
25 uL volume reaction with the following procedures: 
initial denaturation step for 5 min at 95 °C, 35 cycles of 
denaturation for 1 min at 94 °C, annealing for 1 min at 
50 °C for Cyt b and 46°C for COI, extension for 1 min 
at 72 °C. Final extension at 72 °C was conducted for 10 
min. The primers for COI are available from the authors 
upon request. Sequences of Cyt b were amplified with the 
primers 14850-1 (Tanaka-Ueno et al., 1998), CytbB and 
CytbC (Bossuyt and Milinkovitch, 2000). PCR products 
were purified with Gel Extraction Mini Kit, and sequenced 
in both directions with BigDye Terminator Cycle 
Sequencing Kit and ABI PRISM 3730. The sequence data 
were submitted to a BLAST search in GenBank to confirm 
identity. Nucleotide sequences were aligned using Clustal 
X 1.81 (Thompson et al., 1997) with default parameters, 
and then optimized by eye in MEGA 4.0 (Tamura et al., 
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2007). 


2.4 Phylogenetic analyses Genealogical reconstructions 
using Bayesian inference (BI) and maximum parsimony 
(MP) considered COI and Cyt b separately and 
concatenated into a single partition. BI analysis was 
executed in MrBayes 3.1.2 (Ronquist and Huelsenbeck, 
2003) using models estimated with MrModeltest v2.3 
(Nylander, 2004). Consensus frequencies, i. e., Bayesian 
posterior probabilities (BPP), were used to estimate 
nodal support. Four separate runs were performed with 
four Markov chains. Each run was conducted for 3 000 000 
generations and sampled every 100 generations. Log- 
likelihood scores were tracked for stabilization and 
the first 25% of the trees were discarded as burn-in. A 
50% majority rule consensus tree was calculated from 
the remaining trees. MP analyses were implemented 
using PAUP 4.0b10a (Swofford, 2003). The heuristic 
MP searches were executed for 1000 replicates with all 
characters treated as unordered and equally weighted. 
One thousand pseudoreplicates were implemented for 
bootstrap analyses. 


3. Results and Discussion 


3.1 Morphometrics The measurements of all specimens 
were summarized in Table 2. The first three factors 
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explained 87.3% of the total variance. The first factor 
alone explained 81.9% of the variance, with large loading 
for the following 13 characters (the load factor > 0.5): 
SVL, HL, HW, SL, IND, IOD, UEW, FL1, FL2, FL3 FL4, 
FEL, and FTL. The second factor explained another 3.2% 
of the total variance, with the following 15 large loading 
characters: SVL, HL, HW, FEL, IND, IOD, TYE, LAL, 
EL, TL, FTL, FOL, TL3, TL4, and TLS. For the third 
factor, 2.3% variance was explained with large loading 
for the follo-wing 11 characters: SVL, LAL, FL3, FEL, 
FTL, FOL, TL1, TL2, TL3, TL4, and TLS. The plots from 
PC1 and PC2 did not completely define the five species. 
However, the new population from Jiemuxi National 
Nature Reserve (open squares) was clearly resolved and 
well differentiated from R. hanluica (X-pattern; Figure 2). 
This analysis indicated that the two breeding groups did 
not constitute a single species. 


3.2 Genealogy and taxonomy Eighty unique sequences 
were newly determined and deposited in GenBank (Table 
1). The sequences contained about 735 bp of coding 
sequences from Cyt b and 561 bp from COI. No indels 
were observed. Alignment of the concatenated sequences 
contained 1296 bp of which 444 sites were variable and 
385 were potentially phylogenetically informative. 

The MP and BI analyses produced very similar trees; 
only the Bayesian tree was shown (Figure 3). The R. 


[O] R jiemuxiensis 
O R. omeinontis 


x R. hanluica 
-+ R. chaochiaoensis 


* R. zhenhaiensis 


1.00000 0.00000 1.00000 2.00000 3.00000 
PC2 


Figure 2 Principal components analysis (PCA) of 23 morphometric variables for 109 individuals of five species (R. jiemuxiensis sp. nov., R. 
omeimontis, R. hanluica, R. chaochiaoensis, and R. zhenhaiensis). Positive PC 1 scores indicate a proportionately longer trunk, larger head, 
longer snout, larger internarial distance, larger interorbital distance, wider upper eyelid, and longer finger. Positive PC 2 scores correspond to 
a proportionately longer trunk, larger head, longer eye, larger internarial distance, larger interorbital distance, larger tympanum, longer lower- 


arm, longer foot and toe. 
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Figure 3 A Bayesian inference tree based on the combined mitochondrial CO1 and Cyt b nucleotide sequence data. Bayesian posterior 


probabilities > 90% are shown near the respective branches. Localities of corresponding samples are indicated in brackets. * = type locality. 
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longicrus complex, including R. japonica from Japan, was 
resolved as an exclusive lineage with high nodal support. 
Six distinct lineages (I-VI) were recovered within this 
group but with unresolved relationships among them. 
Rana japonica, R. chaochiaoensis, and R. omeimontis 
constituted exclusive lineages I, II, and IV, respectively. 

The two breeding groups from Jiemuxi National 
Nature Reserve were exclusively assigned to lineages 
III and V. There was no mixing of haplotypes within 
the breeding groups at the nature reserve. Lineage V 
exhibited substantial divergence from lineage III and 
the other groups (Figure 2). The null hypothesis of 
conspecificity required non-exclusive maternal lineages 
within the reserve because lineage segregation is not 
expected within a species. The presence of two exclusive 
lineages that breed at different seasons, combined with 
the morphometric analysis, required rejection of the 
null hypothesis and the description of a new species, as 
provided below. 

The genealogical relationships were not only 
informative with respect to the identity of species at 
Jiemuxi National Nature Reserve, but they also had 
additional implications for the validity of some other 
species. Within lineage II], some individuals from Jiemuxi 
National Nature Reserve clustered with R. hanluica and 
R. maoershanensis. Genetic samplings of the named taxa 
are from their respective type localities. The identical 
and nearly identical haplotypes in these three species 
suggested extremely recent isolation events, introgressive 
hybridization, or on-going gene flow. In the first case, 
species recognition was not justified because isolation, if 
present, was likely to be ephemeral; isolation likely owed 
to the current interglacial period. In the hybridization 
scenario, the absence of genetic and morphological 
differentiation did not support this possibility. 
Consequently, our analyses suggested that all members of 
lineage III were best considered to be one species. Based 
on rules of priority in the Code, we recognized lineage III 
as R. hanluica. 

Rana maoershanensis was initially reported as a new 
distribution record of R. chaochiaoensis from Guangxi 
(Lu et al., 2006), but later described as a new species 
(Lu et al., 2007). Fei et al. (2009) did not recognize R. 
maoershanensis and our study suggested it is a junior 
synonym of R. hanluica. This documentation contrasts 
sharply with that of Yang et al. (2011) who used the 
mitochondrial gene 16S to evaluate the status of R. 
maoershanensis. Unfortunately, their analysis recovered 
only one species in the region and the data base contained 
only 58 potentially informative sites for 14 species of 
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frogs. As a consequence of the limited data, most nodes 
at the base of their BI, MP and maximum likelihood 
trees were not supported. This precluded a meaningful 
comparison of our tree with their tree. Nevertheless, their 
report of R. maoershanensis outside of the R. japonica 
complex suggests that the region might contain another 
species of Rana. 

Cryptic species may be common within the complex. 
Rana chaochiaoensis (lineage II) is very distinct (Figure 
2) and this suggests that the species is valid. Although 
an assessment of nuclear genes is required to confirm the 
species’ validity, the association of R. maoershanensis 
with R. chaochiaoensis (Lu et al., 2006) suggests that 
clear morphological diagnoses may not be common in the 
group. The absence of morphological differentiation does 
not preclude recognition of species; cryptic species are a 
reality and their recognition serves both evolutionary and 
conservation biology. 

Our genetic analyses discover that the taxonomy of 
the brown frogs is far more complicated than originally 
imagined. Two genetic lineages that constitute two 
species of brown frogs from Maoershan (III, VI) are 
clearly recovered herein, and these correspond to our field 
observations of two different breeding seasons. Frogs in 
lineage III at Maoershan initiate breeding around October 
although the exact time varies depending on local climate. 
Within this group, breeding can continue to December. We 
have observed tadpoles under the ice during winter, which 
is consistent with the observations of Lu et al. (2007). 
The breeding season of specimens in lineage VI starts in 
February and the same pattern occurs at Yangmingshan. 
Here, it also appears that two cryptic species occur 
sympatrically and this may be a general pattern. 

Lineage VI covers a wide range in China. Specimens 
of R. longicrus, R. zhenhaiensis, and R. culaiensis from 
their respective type localities all nest in this group. Their 
conservative morphology makes the correct determination 
of species challenging, especially in southern China. 
Most often, these species are identified by locality 
data only. Rana zhenhaiensis occurs in most regions 
of southeastern China (Fei et al., 2009). Our analyses 
suggest that specimens of R. zhenhaiensis from Fujian 
and Guangdong provinces might best be assigned to R. 
longicrus because the mtDNA is virtually identical, and 
R. longicrus branches off within R. zhenhaiensis. Further, 
two individuals of R. zhenhaiensis from Hubei and Henan 
cluster with R. culaiensis. Wood frogs from Xinhua, 
Yangmingshan (Hunan), and Maoershan (Guangxi) cluster 
together within lineage VI. Rana zhenhaiensis from 
type locality (Ningbo) forms its own sublineage within 
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lineage VI. Although this pattern is interesting, our sparse 
geographic sampling within lineage VI, and the exclusive 
use of mtDNA sequences, hinders an evaluation of the 
taxonomic status of these species, especially considering 
that sympatric cryptic lineages are common. Previous 
studies based on limited sampling (e.g., Jiang et al., 2001; 
Yang et al., 2001; Che et al., 2007) must be verified 
and collections should be made at different seasons 
of the year. Molecular data should be combined with 
morphology, ecology, and acoustic studies. This approach 
will likely reveal additional cryptic species. 

A more fine-scale sampling from southern China 
should be undertaken in order to understand the diversity 
and relationships among the species of wood frogs in the 
R. longicrus species group. Inclusion of additional genes, 
especially nuclear genes, would be essential. mtDNA data 
alone cannot assess introgression and gene flow. 


4. Taxonomy 


Rana jiemuxiensis sp. nov. (Figures 4, 5) 


Holotype KIZ 05553 (Figure 5), an adult male from 
Jiemuxi National Nature Reserve (28° 52'45” N, 110° 
24' 608” E, elevation 723 m), western Hunan, China. 
Collected by Fang YAN and Zhongmei ZHANG from 10 
to 19 March 2010. 


Paratypes KIZ05245—-KIZ05269, KIZ05270-KI1Z05272, 
KIZ05551, KIZ05554—-K1IZ05556 with the same 
collection date and locality as the holotype. A total of 32 
adult individuals, including 22 males and 10 females. 


Diagnosis Rana jiemuxiensis sp. nov. can be 
distinguished from all other species of the genus Rana by 
the following combination of characters: (1) dorsolateral 
fold from posterior canthus to crotch, slight curving 
above the tympanic membrane; (2) tips of fingers obtuse 
without disks; (3) three metacarpal tubercles in males, 
inner one covered partly by the nuptial pad at the base of 
the first finger, the two outer ones combined at the base 
of the third and forth fingers; (4) gray-whitish nuptial 
pad developed, and divided into three groups; (5) five 
to six wide brownish-black bands at thigh and tibia; (6) 
lineae musculinae absent in males; (7) tibio-tarsal extends 
beyond the snout; (8) breeding season from the end of 
late February to middle March. 


Etymology The species is named for the type locality, 
Jiemuxi, Yuanling County, Hunan. 


Description of holotype SVL 48.48 mm. Head slightly 
longer than broad, snout pointed and projecting; from the 


anterior orbit to snout tip much longer than eye diameter; 
distance between nostril and anterior orbit longer than 
distance between snout-tip and nostril; interorbital 
distance equal to upper eyelid width; tympanum diameter 
slightly shorter than eye diameter; tympanic rim narrow, 
but prominent; loreal region concave, sloping outwards; 
vomerine teeth in short oblique series, anterior edges in 
line with centers of choanae; tongue notched posteriorly, 
not deeply; internal vocal sacs and lineae musculinae 
absent. 

Forearm robust, tips of fingers obtuse, not expanded, 
without circummarginal grooves; relative length of 
fingers: II < IV < I < III; one prominent subarticular 
tubercle on fingers I and I, two small subarticular 
tubercles on fingers III and IV; inner metatarsal tubercle 
strong and larger, ovoid, partly covered by nuptial pad 
at base of finger I; two outer tubercles combined at base 
of fingers III and IV, flat and weak; strong nuptial pad 
divided into three groups, basal one covering anterior 
inner metatarsal tubercle, largest one ventrally extending 
from base of palmar onto medial dorsal surface of finger I 
to anterior edge of subarticular tubercle, the smallest one 
extending from anterior of largest nearly to tip of finger; 
nuptial pad with off-white dense spinules. 

Hindlimb long, tibio-tarsal joint beyond snout-tip; 
indistinct skin ridge on dorsum of tibia; heels overlapping 
more when legs are held perpendicular to body; tips of 
toes similar to fingers; relative length of toes: I < H <HI< 
V < IV; toes webbed fully, the web of toe IV reaching 
the second joint from tip; subarticular tubercles distinct, 
especially joint of toe and metatarsus; inner metatarsal 
tubercle andovoid, small but well developed; outer 
metatarsal tubercle indistinct. 

Skin rather smooth; some small granules on body 
side or near vent; temporal fold distinct, a somewhat 
triangular gray patch behind eye and anterior to temporal 
fold; thickened dorsolateral fold from posterior canthus 
to crotch, slightly curved at upper tympanum and 
crossing temporal fold; elongate rectal gland reaching 
temporal fold and expanding; ventral surface smooth; 
some somewhat whitish granules on the posterior ventral 
surface of femurs. 


Coloration Dorsum of body and body side grayish brown 
(gray in ethanol), with scattered black spots; the side of 
snout and upper lip before tympanum brown (grayish 
black in ethanol); lower lip whitish with black spots; 
forelimb grayish brown (gray in ethanol), two indistinct 
black bars on dorsum of each arm; dorsum of hand with 
white areas; dorsum of thigh and tibia grayish brown (gray 
in ethanol), with four and three black bars, respectively; 
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Figure 4 (a) Lateral view of mating R. jiemuxiensis on land during the breeding season. (b) Lateral view of mating R. jiemuxiensis under 
water during the breeding season. (c) Lateral view of a male R. jiemuxiensis. (d) Lateral view of a female R. jiemuxiensis. 


Figure 5 Holotype specimen of R. jiemuxiensis. (a) Dorsal view. (b) Ventral view. (c) Details of right hand showing the nuptial pad. (d) 
Details of right foot showing the extent of webbing. 
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Figure 6 Habitat of the type locality of R. jiemuxiensis (a); 
clutches of eggs (b, c). 


anterior thigh and leg with black spots; tarsus and foot 
side grayish brown (gray in ethanol), tarsus with three 
black bars; dorsum of foot whitish with black mottling; 
throat and belly whitish with black mottling. 


Habitat and breeding behaviour The breeding season 
is from February to March with a peak in early March. 
The exact timing can be influenced by either precipitation 
or human agricultural activities. Breeding individuals 
grouped at a nearby rice paddy, and the water came 
from either rainfall or irrigation (Figure 6). Since 2008, 
breeding lasted only for about ten days. The number 
of attending males always exceeded that of females. 
Males, which possess an obvious nuptial pad, always 
occupied a location to attract females with calls. Egg 
masses contained 200-1500 eggs. No parental care was 
observed. In 2010, froglets appeared in May, suggesting 
that development to metamorphosis required about two 
months. 
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